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^ ■ 1 Introduction 



Recently an increasing interest has been attracted by the topic of rotationally real invariant ran- 
dom fields on the sphere due to applications to the statistical analysis of Cosmological and 
Astrophysical data (see [MP04], [Mar06] and [AK05]). 

Some results concerning their structure and their spectral decomposition have been obtained in 
[BM07], where a peculiar feature has been pointed out, namely that if the development into spherical 
harmonics 

oo m 

T = ^ ^ at,mXt, m 
e=i -m 

of a rotationally invariant random field T is such that the coefficients a^ m , £ = 1,2, . . . ,0 < m < £ 
are independent, then the field is necessarily Gaussian (the other coefficients are constrained by the 
condition a^_ m = (— l) m a^ m ). This fact (independence of the coe:fficients+isotropy=>Gaussianity) 
is not true for isotropic random fields on other structures, as the torus or Z (which are situations 
on which the action is Abelian). 

This property implies in particular that non Gaussian rotationally invariant random fields on the 
sphere cannot be simulated using independent coefficients. 
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In this note we show that this is a typical phenomenon for homogeneous spaces of compact non- 
Abelian groups. This should be intended as a contribution to a much more complicated issue, i.e. 
the characterization of the isotropy of a random field in terms of its random Fourier expansion. 

In §2 and 3 we review some background material on harmonic analysis and spectral repre- 
sentations for random fields. §4 contains the main results, whereas we moved to §5 an auxiliary 
proposition. 

2 The Peter- Weyl decomposition 

Let 1" be a compact topological space and G a compact group acting on 2£ transitively. We 
denote by mc the Haar measure of G. We know that there exists on if a probability measure m 
that is invariant by the action of G, noted x — ► g~ x x, g e G. We assume that both m and mo arc 
normalized and have total mass equal to 1. We shall write L 2 (%~) or simply L 2 instead of L 2 ( ,ra). 
Unless otherwise stated the spaces L 2 are spaces of complex valued square integrable functions. We 
denote by L g the action of G on L 2 , that is L g f(x) = f(g~ 1 x). 

Let 3£ be the set of equivalence classes of irreducible unitary representations of G which occur 
in the decomposition of L 2 (i?T, m). Since the action of G commutes with the complex conjugation 
on L 2 (i?T,m), it is clear that for any irreducible subspace H, H, its conjugate subspace is also 
irreducible. If H = H, we can find orthonormal bases (<j>k) for H which are stable under conjugation; 
for instance we can choose the <j>k to be real. If H ^ H, then there are two cases according as the 
action of G on H is, or is not, equivalent to the action on H . If the two actions are inequivalent, then 
automatically H _L H. If the actions are equivalent, it is possible that H and H are not orthogonal 
to each other. In this case H D H = as both are irreducible and S = H + H is stable under G and 
conjugation. In this case we can find K C S stable under G and irreducible such that K _L K and 
S = K ® K is an orthogonal direct sum. The proof of this is postponed to the Appendix so as not 
to interrupt the main flow of the argument. We thus obtain the following orthogonal decomposition 
of L 2 (3£ ,m), compatible with complex conjugation: 

L 2 (,r,m) = H t ® {H t ®Hl) (2.1) 

where the direct sums are orthogonal and 

i e J° ^> Hi—Hi, ie Hi LHi. 

We can therefore choose an orthonormal basis (4>ik) for L 2 (££ , m) such that for i <G J^°, (4>ik)i<k<di 
is an orthonormal basis of Hi stable under conjugation, while, for i G (4>ik)i<k<di is & n or- 

thonormal basis for Hi, where di is the dimension of Hf, then, for i G {4>ik)i<k<di is an 

orthonormal basis for Hi. Such a orthonormal basis {<f)ik)ik of L 2 (%~,m) is said to be compatible 
with complex conjugation. 

Example 2.1 S£ = S 1 , the one dimensional torus. Here G = Z and Hk, k € Z is generated by the 
function 7^(0) = e lk0 . Hk = H_k an d Hk -L Hk for k ^ 0. All of the H^s are one-dimensional. 

Recall that the irreducible representations of a compact topological group G are all one-dimensional 
if and only if G is Abelian. 

Example 2.2 G = 5*0(3), ^ = S 2 , the sphere. A popular choice of a basis of L 2 (^T,m) are the 
spherical harmonics, (l£, m )-^< m <£, I <= N (see [VK91]). Hi — span((Y^ m )_£< m <^) are subspaces of 
L 2 (%~,m) on which G acts irreducibly. We have Ye, m — (— l) m y«,-m and Y^o is real. 

By choosing <j>e lTn — Ye. m for m > and 4>e im = {—^-) m Yl,m for m < 0, we find a basis of Hi such 
that if 4> is an element of the basis, then the same is true for <j>. Here dim(Hi) = 21+1, Hi = Hi, 
so that in the decomposition (2.1) there arc no subspaces of the form Hi for J* + . 
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3 The Karhunen-Loeve expansion 

We consider on X a real centered square integrable random field (T(x)) X £x ■ We assume that there 
exists a probability space (£l,&,P) on which the r.v.'s T(x) are defined and that (x,co) — > T(x, w) 
is ^ (,$£") (g) measurable, ^"(^T) denoting the Borel cr-field of SC . We assume that 



E 



/ T{xf dm(x)} = M < +00 (3.2) 
J X i 



which in particular entails that x — > T x (w) belongs to L 2 (m) a.s. Let us recall the main elementary 
facts concerning the Karhuncn-Locvc expansion for such fields. We can associate to T the bilinear 
form on L 2 (m) 

T(f,g) = E\[ T(x)f(x)dm(x) f T(y)g(y) dm(y)} (3.3) 
J X J X J 

By (3.2) and the Schwartz inequality one gets easily that 

\T(f,g)\<M\\f\\ 2 \\g\\ 2 . 

Therefore, by the Riesz representation theorem there exists a function R € L 2 (^ x % ', m® m) such 
that 

T(f,g)= / f(x)g(y)R(x,y)dm(x)dm(y) . 
We can therefore define a continuous linear operator R : L 2 (m) — > L 2 (m) 



Rf(x)= / R{x,y)f{y)dm{y) . 

Jar 

It can be even be proved that the linear operator i? is of trace class and therefore compact (see 
[Par05] for details). Since it is self-adjoint there exists an orthonormal basis of L 2 (^,m) that is 
formed by eigenvectors of R. 

Let us define, for <p <G L 2 { 3£ , m), 



a(</>) = / T(x)(/>(x) rfm(i) , 



Let A be an eigenvalue of R and denote by E\ the corresponding eigenspace. Then the following is 
well-known. 

Proposition 3.3 Let <fr E E\. 

a) If ip <E L 2 (%~,m) is orthogonal to (f>, a(ip) is orthogonal to a{(f>) in L 2 (f2,P). Moreover 

E[Kvoi 2 ] = A||v>||i. 

b) If <ft is orthogonal to <\>, then the r.v.'s 3?a(0) and 9a(0) are orthogonal and have the same 
variance. 

c) If the held T is Gaussian, a((j>) is a Gaussian r.v. If moreover <j> is orthogonal to 4>, then a((f>) 
is a complex centered Gaussian r.v. (that is 3?ai and Soi are centered, Gaussian, independent and 
have the same variance). 

Proof, a) We have 

E[a(cf))a(ip)} = E[ / T(x)<t>(x)dm{x) ( T(y)^(y) dm(y) = 

J X J X J 

= / R(x,y)<j>(x)ip(y)dm(x)dm(y) = X <j>(y)ip(y)dm(y) = \(<j>,ip) . 
JXxX J X 
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From this relation, by choosing first tp orthogonal to (f> and then ip — (p, the statement follows. 

b) From the computation in a), as a(<j>) — a{(j>), one gets E[a(</>) 2 ] = \{4>,<j>). Therefore, if <j> is 
orthogonal to <fi, E[a(0) 2 ] = which is equivalent to $la(<fr) and 9a(</>) being orthogonal and having 
the same variance. 

c) It is immediate that a{<f>) is Gaussian. If 4> is orthogonal to <f>, a(<p) is a complex centered 
Gaussian r.v., thanks to b). 

■ 

If {<t>k)k is an orthonormal basis that is formed by eigenvectors of R, then under the assumption 
(3.2) it is well-known that the following expansion holds 

oo 

T(x) = J2 a (MMx) (3.4) 
fe=i 

which is called the Karhunen-Loeve expansion. This is intended in the sense of L 2 ($y , to) a.s. in uj. 
Stronger assumptions (continuity in square mean of x — > T(x), e.g.) ensure also that the convergence 
takes place in L 2 {£1, P) for every x (see [SW86], p. 210 e.g.) 

More relevant properties are true if we assume in addition that the random field is invariant 
by the action G. Recall that the field T is said to be (weakly) invariant by the action of G if, for 
fi, ■ ■ ■ fm € L 2 (^) the joint laws of (T(/i), . . . , T(.f m )) and (T(L s /i), . . . , T((L g f m )) are equal for 
every g € G. Here we write 

T(f) = f T{x)f(x) dm(x), f e L 2 (^) . 

If, in addition, the field is assumed to be continuous in sqare mean, this imples that for every 
x\, . . . ,x m <G (T(xi), . . . , T(x m )) and (T(g x\), . . . 1 T(g~ 1 x m )) 1 have the same joint laws for 
every g e G. If the field is invariant then it is immediate that the covariance function R enjoys the 
invariance property 

R(x,y) = R(g~ 1 x,g~ 1 y) a.e. for every g e G (3.5) 

which also reads as 

L g (Rf) = R(L g f) . (3.6) 

Then, thanks to (3.6), it is clear that G acts on E\. Therefore E\ is the direct sum of some of 
the Hi's introduced in the previous section. Moreover it is a finite direct sum, unless A = 0, as the 
eigenvalues of a compact operator that are different from cannot have but a finite multiplicity. It 
turns out therefore that the basis {4*ik)ik of L 2 introduced in the previous section is always formed 
by eigenvectors of R. 

Moreover, if some of the H^s are of dimension > 1, some of the eigenvalues of R have necessarily a 
multiplicity that is strictly larger than 1. As pointed out in §2, this phenomenon is related to the non 
commutativity of G. For more details on the Karhunen-Loeve expansion and group representations 
see [PP05]. 

Remark that if the random field is isotropic and satisfies (3.2), then (3.4) follows by the Peter- 
Weyl theorem. Actually (3.2) entails that, for almost every u>, x — > T(x) belongs to L 2 (^~,m). 

Remark 3.4 An important issue when dealing with isotropic random fields is simulation. In this 
regard, a natural starting point is the Karhunen-Loeve expansion: one can actually sample random 
r.v.'s a((f>k), (centered and standardized) and write 

n 

T n{x) = ^2 V^fc a(<j> k )4>k (3.7) 
fc=i 
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where the sequence (Xk)k is summable. The point of course is what conditions, in addition to those 
already pointed out, should be imposed in order that (3.7) defines an isotropic field. In order to 
have a real field, it will be necessary that 

a{4> k ) = a^u) (3.8) 

Our main result (see next section) is that if the a((/>fc)'s are independent r.v.'s (abiding nonetheless 
to condition (3.8)), then the coefficients, and therefore the field itself are Gaussian. 

If Hi C L 2 ( JT, m) is a subspacc on which G acts irreducibly, then one can consider the random field 

Iff, (*) (*) 

where the <pj are an orthonormal basis of Hi. As remarked before, all functions in Hi are eigenvectors 
of R associated to the same eigenvalue A. 

Putting together this fact with (3.4) and (2.1) we obtain the decomposition 

Example 3.5 Let T be a centered random field satisfying assumption (3.2) over the torus T, whose 
Karhuncn-Locvc expansion is 

+ oo 

T(8) = Y j a k e lke 1 9 eT . 

— oo 

Then, if T is invariant by the action of T itself, the fields (T(8))g and (T(9+9'))g are equi-distributed, 
which implies that the two sequences of r.v.'s 

(afc)-oo<fc<+oo and (e tke a k )- oc<k<+QO (3.9) 

have the same finite distribution for every 9' e T. Actually one can restrict the attention to the 
coefficients (afc)o<fc<+oo, as necessarily o_fe = a k - 

Conversely it is clear that if the two sequences in (3.9) have the same distribution for every 
9' e T, then the field is invariant. 

Condition (3.9) implies in particular that, for every k, — oo < k < +oo, k ^ the distribution of 
a k must be invariant by rotation (i.e. by the multiplication of a complex number of modulus 1). 

If one assumes moreover that the r.v.'s a k , are independent, then every choice of a distribution 
for a/j.,0 < k < +oo that is rotationally invariant gives rise to a random field that is invariant with 
respect to the action of T. 

4 Independent coefficients and non-Abelian groups 

In this section we prove our main results showing that, if the group G is non commutative and 
under some mild additional assumptions, independence of the coefficients of the Fourier development 
implies their Gaussianity and, therefore, also that the random field must be Gaussian. We stress 
that we do not assume independence of the real and imaginary parts of the random coefficients. 

Proposition 4.6 Let !% be an homogeneous space of the compact group G. Let Hf C L 2 (%~,m) 
be a subspace on which G acts irreducibly, having a dimension > 2 and such that if f G H^~ then 
f H^~ . Let {<pk)k be an orthonormal basis of H^ and consider the random held 

T h+( x ) =E a A( i ) • 

k 
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for a family of r.v.'s (cik)k C L 2 (fl, P). Then, if the r.v.'s cii are independent, the field T H + is 
G-invariant if and only if the r.v.'s (dk)k are jointly Gaussian and E(|afc | 2 ) = c (and therefore also 
the held T H + is Gaussian). 



Proof. Since G acts irreducibly on 



we have 



4>k{g 1 x) = ^D M (g)<^(x) , 




T{g- 1 x)=Y,at<j>i{x) 



where 



di 

at = D k j(g)a k . 
k=i 



If the field is G-invariant, then the vectors {at)t have the same joint distribution as (a k )k and in 
particular the (at)t are independent. One can then apply the Skitovich-Darmois theorem below (see 
[KLR73] e.g.) as soon as it is proved that g £ G can be chosen so that D k j(g) ^ for every k,£. 
This will follow from the considerations below, where it is proved that the set Z k ,i of the zeros of 
Dk t e has measure zero. 

Indeed, let G\ be the image of G in the representation space so that G\ is also a connected 
compact group, and is moreover a Lie group since it is a closed subgroup of the unitary group \J(di). 
If the representation is non trivial, then G\ ^ {1} and in fact has positive dimension, and the D k ,i 
are really functions on G\. For any fixed k, t the irreducibility of the action of G\ implies that D k ,t 
is not identically zero. Indeed, if this were not the case, we must have {g4>t, 4>k) = for all g £ G\, 
so that the span of the g4>i is orthogonal to (f> k ; this span, being Gi-invariant and nonzero, must be 
the whole space by the irreducibility, and so we have a contradiction. 

Since D k e is a non zero analytic function on G\, it follows from standard results that Z k t has 
measure zero. Hence lj fc ^ Z k i has measure zero also, and so its complement in G\ is non empty. 

■ 

We use the following version of the Skitovich-Darmois theorem, which was actually proved by 
S. G. Ghurye and I. Olkin [G062] (see also [KLR73]). 

Theorem 4.7 Let X\, . . . , X r be mutually independent random vectors in R n . If the linear statistics 



are independent for some real nonsingular n x n matrices Aj,Bj, j = 1, . . . ,r, then each of the 
vectors X\, . . . ,X r is normally distributed. 

We now investigate the case of the random field Th, when H is a subspace such that H = H. In 
this case we can consider a basis of the form 4>-k, ■ ■ ■ , 4>k, k < I, with <j)-k — 4>k- The basis may 
contain a real function </> , if dim if is odd. Let us assume that the random coefficients a k , k > 
are independent. Recall that a- k — TTj:. 



r 



r 



Li - y^AjXj, 



L 2 - ^2 B i X i 



G 



The argument can be implemented along the same lines as in Proposition 4.6. More precisely, if 
mi > 0, m 2 > 0, the two complex r.v.'s 

£ 

m =~ e (4.10) 

m——t 

have the same joint distribution as a mi and a„ l2 . Therefore, if mi 7^ m2, they arc independent. 
Moreover a_ m = a^, so that the previous relation can be written 

i 

a mi = D Qmi (g)a + [D m , mi (g)a m + D_ m ^ ni (g)a^j 

m—l 

t 

a m2 = D , m2 {g)a + (D m . m2 (g)a m + D- m . m2 (g)a^j 

m—l 

In order to apply the Skitovich-Darmois theorem, we must ensure that g G G can be chosen so that 
the real linear applications 

z -> D miTni (g)z + D_ mimi (g)z, m = 1, ...,£, i = 1, 2 (4.11) 

are all non singular. It is immediate that this condition is equivalent to imposing that \D m _ mi (g) \ =fi 
\D- m , m( (g)\. 

We show below that (4.11) is satisfied for some well-known examples of groups and homogeneous 
spaces. We do not know whether (4.11) is always satisfied for every compact group. We are therefore 
stating our result conditional upon (4.11) being fulfilled. 

Assumption 4.8 There exist g £ G, < mi < m 2 < £ such that 

—m.rrii (g)\ 

for every < m < £. 

We have therefore proved the following. 

Proposition 4.9 Let X be an homogeneous space of the compact group G. Let Hi C L 2 (^",m) 
be a subspace on which G acts irreducibly, having a dimension d > 2 and such that Hi = Hi. Let 
(4>k)k be an orthonormal basis of Hi such that <j)~k — 4>k an d consider the random field 

T Hi(x) =^a k (j) k (x) 

k 

where the r.v.'s a kl k > are centered, square integrable, independent and a_fc = a^. Then Tf^ is 
G-invariant if and only if the r.v. 's {a k )k>o are jointly Gaussian and E(\ak\ 2 ) = c (and therefore also 
the field THi is Gaussian). 

Putting together Propositions 4.6 and 4.9 we obtain our main result. 

Theorem 4.10 Let 0£ be an homogeneous space of the compact group G. Consider the decompo- 
sition (2.1) and let ((4>ik)iej?° ,{(f>ik,4'ik)ieJ^+) be a basis of L 2 (G) adapted to that decomposition. 
Let 

T= ^ ^^ a ik4>ik+ ^ ^ {aik4>ik + aik<f>ik) 

k i£j>+ k 
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be a random field on S£ , where the series above are intended to be converging in square mean. Assume 
that T is isotropic with respect to the action of G and that the coefficients (a i/ t) ie ,^o fe > , (ajfc)j £ /+ 
are independent. If moreover 

a) the only one- dimensional irreducible representation appearing in (2.1) are the constants; 

b) there are no 2 -dimensional subspaces H C L 2 (5£), invariant under the action of G and such 
that ~H = H; 

c) The random coefficient corresponding to the trivial representation vanishes. 

d) For every H C L 2 (%~), irreducible under the action of G and such that H — H , Assumption 
4-8 holds. 

Then the coefficients (aik)iey° ,fc>(b ( a ik)ie.j?+ are Gaussian and the field itself is Gaussian. 

Let us stress with the following statements the meaning of assumption a)-d). The following result 
gives a condition ensuring that assumption b) of Theorem 4.10 is satisfied. 

Proposition 4.11 Let U be an irreducible unitary 2-dimensional representation of G and let Hi 

and H 2 be the two corresponding subspaces of L 2 (G) in the Peter- Weyl decomposition. Then if U 
has values in SU{2), then Hi = H 2 ^ Hi. 

Proof. If we note 




then one can assume that Hi is generated by the functions a and c, whereas H 2 is generated by b 
and d. It suffices now to show that a is orthogonal both to a and c. But, the matrix U(g) belonging 
to SU(2), we have a{g) = d{g) E H 2 . 

■ 

Recall that the commutator Go, of a topological group G is the closed group that is generated by 
the elements of the form xyx~ l y~ l 

Corollary 4.12 Let G be a compact group such that its commutator Go coincides with G himself. 
Then assumptions a) and b) of Theorem 4-10 are satisfied. In particular these assumptions are 
satisfied if G is a semisimple Lie group. 

Proof. Recall that if Go = G, G cannot have a quotient that is an abelian group. If there was a 
unitary representation U with a determinant not identically equal to 1, then g — > det(£7(G)) would 
be an homomorphism onto the torus T and therefore G would possess T as a quotient. The same 
argument proves that G cannot have a one dimensional unitary representation other than the trivial 
one. One can therefore apply Proposition 4.11 and b) is satisfied. 

■ 

Remark 4.13 It is easy to prove that Assumption 4.8 is satisfied when 3£ = §> 2 and G = 50(3). 
As mentioned in [BM07], this can be established using explicit expressions of the representation 
coefficients as provided e.g. in [VMK88]. 

In the same line of arguments it is also easy to check the same in the cases X — 50(3), 
G = 50(3) and JT = SU(2), G = SU(2). 

Remark 4.14 As far as condition c) of Theorem 4.10, let us remark that the coefficient of the trivial 
representation corresponds to the empirical mean of the field. As any random field can be decom- 
posed into the sum of its empirical mean plus a field whose coefficient of the trivial representation 
vanishes, our result can be interpreted in terms of Gaussianity of this second component. 
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5 Appendix 



Proposition 5.15 Let V be a finite dimensional Hilbert space on which G acts unitarily, and let V 
be equipped with a conjugation a(v — > v) commuting with the action of G. Let H be an irreducible 
G-invariant subspace and let V = H + H . 

a) If the actions of G on H and H are inequivalent, then H _L H and V — H © H . 

b) If the actions of G on H and H are equivalent, then either H = H or we can find an 
irreducible G-invariant subspace KofV such that K _L K and V — K © K . 

Proof. Let P be the orthogonal projection V — > H and A its restriction to H. Then, for every h € H, 
hi e H and g e G, we have 

(g(Ah),ti) = (Ah,gh') = (h,gti) = (gh,ti) = (A(gh),h') 

From this we get that G acts on A(H). The action of G on H being irreducible, we have either 
A(H) = {0} or A(H) = H. In the first case H is already orthogonal to H . Otherwise A intertwines 
the actions on H and on H, so that these are equivalent and V — H © H . 

V being the sum of two copies of the representation on H, there is a unitary isomorphism 
V ~ H © C 2 where C 2 is given the standard scalar product. So we assume that V = H © C 2 . G 
acts only on the first component, so that G acts irreducibly on every subspace of the form H © Z , 
Z being a one dimensional subspace of C 2 . 

Let us identify the action of a on H ©C 2 . Let cro be the conjugation on V defined by <7o (it ©v) = 
u© v where v —> v is the standard conjugation (z\,z-i) — > (zT, zj). Then ao-Q is a linear operator 
commuting with G and so is of the form 1 © L where L(C 2 — > C 2 ) is a linear operator. Hence 

a(u © v) = aero (w ®v) = u © Lv. 

If Z is any one dimensional subspace of C 2 , H © Z is G-invariant and irreducible, and we want to 
show that for some Z,H®Z±H®Z a , i.e., Z _L Z a . Here Z a = a(Z). 

For any such Z, let v be a nonzero vector in it; then the condition Z _L Z a becomes (v, Lv) = 
where (, ) is the scalar product in C 2 . Since (, ) is Hermitian, B(v, w) :— (v, Lw) is bilinear and we 
want v to satisfy B(v, v) — 0. This is actually standard: indeed, replacing B by B + B T (which just 
doubles the quadratic form) we may assume that B is symmetric. 

If B is degenerate, there is a nonzero v such that B(v, w) = for all w, hence B(v, v) = 0. If B 
is nondegenerate, there is a basis vi,v 2 for C 2 such that B(vi,vj) = Sij. Then, if w = v\ +iv 2 where 
i = V^T, B(w, w) = 0. 
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